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The rate of second layer nucleation - the formation of a stable nucleus on top of a two-dimensional 
island - determines both the conditions for layer-by-layer growth, and the size of the top terrace 
of multilayer mounds in three-dimensional homoepitaxial growth. It was recently shown that con- 
ventional mean field nucleation theory overestimates the rate of second layer nucleation by a factor 
that is proportional to the number of times a given site is visited by an adatom during its residence 
time on the island. In the presence of strong step edge barriers this factor can be large, leading to a 
substantial error in previous attempts to experimentally determine barrier energies from the onset 
of second layer nucleation. In the first part of the paper simple analytic estimates of second layer 
nucleation rates based on a comparison of the relevant time scales will be reviewed. In the main part 
the theory of second layer nucleation is applied to the growth of multilayer mounds in the presence of 
strong but finite step edge barriers. The shape of the mounds is obtained by numerical integration 
of the deterministic evolution of island boundaries, supplemented by a rule for nucleation in the 
top layer. For thick films the shape converges to a simple scaling solution. The scaling function is 
parametrized by the coverage 9 C of the top layer, and ta kes the fo rm of an inverse error function cut 
off at 6 C . The surface width of a film of thickness d is y/ (1 — 9 c )d. Finally, we show that the scaling 
solution can be derived also from a continuum growth equation. 



I. INTRODUCTION 



An atom adsorbed on a crystal surface can move from one atomic layer to another (usually from a higher layer to 
a lower one) only by crossing a step. Irrespective of whether it occurs by simple hopping or by concerted exchange, 
step crossing is often associated with an additional energy barrier AEg , which reduces the effective rate of interlayer 
transport, v' , below the rate v for adatom diffusion within one atomic layer jj],^]]. Here we are assuming conventional 
Arrhcnius laws for v and v' ', 

V = V e-^/teT v l = v , q e -E s /k B T (1) 

and define AEs — Eg — Ed- The consequences for multilayer growth are twofold. On the one hand, the confinement 
of atoms on top of the two-dimensional islands which form during the growth of the first layer increases the nucleation 
rate of the second layer, which therefore may take place before the first layer islands have coalesced, leading to three- 
dimensional growth H . On the other hand, once three-dimensional mound- like features have formed through this (or 
some other) mechanism, the preferential attachment of atoms to the ascending step of the vicinal terraces on the sides 
of the mounds implies an effective uphill mass current, which stabilizes and amplifies the mound formation (see 



Section IV). The ubiquity of step edge barriers on metal surfaces implies that also three-dimensional mound growth 
should be ubiquitous, and indeed the phenomenon has been observed, among other systems, in the homoepitaxy of 
Pe(100) |,|, Rh(lll) §, Cu(100) || and Pt(lll) @. 

These patterns evoke questions about the kinetic processes that determine the shape of individual mounds |ll| , |l2| , 
as well as about those processes that contribute to the coarsening of the pattern by transferring mass from smaller 
mounds to larger ones |13|-[l6|j. In the present contribution we focus on the first set of questions. Starting from recent 
insights into the process of nucleation on top of islands p7|-p0[, we develop a minimalistic model for the growth of 
a single mound, which allows us to analytically determine the asymptotic mound shape, as well as to corroborate 
earlier results that link the size of the topmost terrace of the mound to the rate of interlayer transport [|l7j . In the 
last section, the asymptotic mound shape is obtained as the solution of a continuum equation for the growing surface. 
This extends a previously established connection between macroscopic and atomistic levels of description [[ll) to a 
somewhat more realistic situation. 



II. THE RATE OF SECOND LAYER NUCLEATION 
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A. STRONG BARRIERS, IRREVERSIBLE AGGREGATION 
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FIG. 1. Kinetic processes involved in second layer nucleation. The figure shows a circular island viewed from the side. 

We consider the geometry illustrated in Figure 2. Atoms are being deposited at rate F onto a circular island of 
radius R. They diffuse on the island with an in-layer hopping rate v, and descend from it with an (average) interlayer 
hopping rate v 1 < v. We are looking for the probability per unit time, u>, for a nucleation event to occur on the island. 
Assuming that dimers of adsorbed atoms are stable (as they will be at sufficiently low temperatures), nucleation takes 
place as soon as two adatoms meet. 

The kinetic rates F, v and v 1 combine to form three relevant time scales^] The interarrival time 

At = —^r- (2) 
ttR 2 F y ' 

between subsequent depositions onto the island; the diffusion time 

t d ~ R 2 /v (3) 

required for an atom to diffuse once across the island, or for two atoms to meet; and the residence time r which an 
atom spends on the island before descending, when no nucleation event takes place. The calculation of the residence 
time requires the solution of a stationary diffusion problem with appropriate boundary conditions at the step edge, 
which yields |Q 



T = 



^« 2 + f)' <«> 



The residence time is comparable to the diffusion time td if the suppression of interlayer transport is weak, in the 
sense that v jv' <C R, while in the opposite regime of strong step edge barriers, v jv' 3> R, the residence time becomes 
t = R/2v' independent of v. 

We will focus on the latter regime in the following. Then r 3> rrj, which implies that two atoms are certain to 
meet once they are present simultaneously on the island. The probability for this to occur is r/(r + At) w r/At if 
At 3> t, which is true for reasonable deposition fluxes. Multiplying this with the total number of atoms deposited 
onto the island per unit time, we obtain the expression JTtJ 

T 7T 2 F 2 i? 5 

W = jMf = ^ r ~ (5) 
for the nucleation rate, which is exact under the stated conditions. 

B. COMPARISON WITH MEAN FIELD THEORY AND EXPERIMENTS 
Equation (||) does not agree with the expression 



1 Throughout lengths will be measured in units of the distance between adsorption sites on the surface. 
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7T aF 2 iyR i 
f ~ 16 



(6) 



derived from conventional (mean field) nucleation theory ||] ; here a denotes a dimensionless capture number of order 
unity. In fact u> m { exceeds the true nucleation rate by a factor which is of the order of (vjv')jR ^> 1. The reason 
can be traced back to the fact that mean field nucleation theory treats the diffusing atoms as noninteracting, which 
implies that they continue to diffuse even after having formed a dimer, thus accumulating too many nucleation events; 
in general the ratio lo^/lo is proportional to the mean number of times a site on the island is visited by an adatom 
during its lifetime p0| . 

In the analysis of experimental data, where in the past the expression (^) has mostly been used instead of (|^), 
this overcounting implies that the value of the step edge barrier AE§ tends to be underestimated. For example, for 
Pt/Pt(lll) with CO-decorated steps, the value AE§ ~ 0.36 eV was obtained using (§), while the mean field expression 
(0) yields only 0.28 eV p7|| . A reanalysis of a second layer nucleation study on Ag(lll) gave AEs ~ 0.15 eV 
instead of the estimate 0.12 eV reported in the original publication, provided the preexponential factor in (|l|) was 
assumed to have the conventional value v' = 10 13 s _1 . The attempt to determine both the step edge barrier and the 
prefactor by fitting the Arrhenius form ([!]) to data obtained at two different temperatures produces unphysically large 
values for z/q, which indicates that the data are not sufficiently accurate for this purpose p2fl. 



C. REVERSIBLE AGGREGATION 



To some extent the above derivation of the nucleation rate can be extended to the case of reversible aggregation, 
where a minimum number of i* + 1 > 2 atoms have to come together to form a stable nucleus Jl8| , ^9[ . The encounter 
time required for the atoms to meet is then no longer of the order of the diffusion time td, but is (under the assumption 
that no metastable clusters exist) given by 

T cnc ~R 2l '/v. (7) 

For i* > 1 it is therefore possible that td <C t <C t ci1c , which implies that nucleation is limited by the encounter of 
the atoms, rather than by their simultaneous presence on the island. The probability p cnc that the i* + 1 atoms, once 
on the island, meet before one of them escapes is then of the order of r /r cnc <C 1, while p cnc ~ 1 when r 3> r cnc . 
In general, the nucleation probability for a freshly deposited atom is proportional to the product of p cnc and the 
probability that i* atoms are already present on the island. Since the latter is of the order of (t/ At) 1 , the nucleation 
rate can be written as 

Pcnc- (8) 



(At) < *+ 1 



Interestingly, in the encounter limited regime where p enc <C 1, this is found to coincide in order of magnitude with 
the prediction of mean field theory. The evaluation of (^) in the two regimes defined by the possible orderings of the 
relevant time scales (to <C T cnc < t < At or < t < T cnc <§; At) then yields the expressions 



where 




Equations (||) and ( |l0| ) respectively describe the scaling regimes II and III of In addition, a nonstationary 

regime I exists, where the residence time is effectively infinite and no stationary balance between deposition and 
escape of atoms from the island is established, as well as a regime IV corresponding to weak step edge barriers. These 
regimes play no role in the context of the present paper. 
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III. THE WEDDING CAKE MODEL 



Our model for the shape of a single mound is a stack of concentric, circular islands (Figure 3). The base is an island 
of radius Rq which does not grow, and the radius of the n'th island is denoted by R n . 



n=n top 



n=l 



*0 



FIG. 2. A stack of concentric circular islands seen from the side. 

A. RATE EQUATIONS FOR LAYER COVERAGES 

The rate at which the n'th island in the stack grows can be computed from the solution of the stationary diffusion 



equation for the annular region separating two islands |23|. We simplify the problem even further, and set the 



interlayer transport rate on the sides of the mound to zero. This is motivated by the fact (to be demonstrated below) 
that the mound steepens indefinitely, and hence the annular regions separating two islands become very narrow. The 
probability for an atom deposited in such an annular region of width I to attach to the ascending, rather than the 
descending step is given by fl24j ] 

l/2 + v/v'l , s 

1 + v/v'l 

which tends to unity for I <C v /v' . Thus interlayer transport will be increasingly suppressed as the mound steepens, 
and it seems reasonable to neglect it from the outset^]. 

In the absence of interlayer transport the layer coverages 6 n — (R n /Ro) 2 satisfy the simple, linear equations [jll 26 



^ = F(9 n . 1 -9 n ), (13) 
at 

which express the fact that all atoms landing on top of the exposed part of layer n — 1 attach to the boundary of layer 
n. For the top layer n top descent of atoms cannot be neglected - since there are no sinks for atoms on the top layer, 
all atoms deposited there have to attach to the descending step. The top terrace therefore absorbs all atoms landing 
on layers n top and n top — 1, and grows according to 

dOn, 

—H^L = F9 nt 14 
dt ,op V ' 



B. NUCLEATION RULES 

Equations (Jl^JlJ) have to be supplemented by a rule for the nucleation of a new top terrace. A simple choice would 
be to posit that nucleation occurs whenever the current top layer reaches some critical coverage p3[ | . More realistically, 
nucleation should be treated as a stochastic process governed by the nucleation rate to computed in Section |J. A 
deterministic rule which is close in spirit to stochastic nucleation can be obtained as follows. Suppose the current top 
terrace has nucleated at time t = 0, and denote its radius by Rtop{t)- Then the probability that no nucleation has 
occurred on the top terrace up to time t is given by (l7[ 



2 In addition, setting u' — makes it impossible for the boundaries of different islands in the stack to collide, which otherwise 
cannot generally be ruled out |25|. 
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Po(t) = cxp (- J dt' uj(R top (t'))^j , (15) 

and the probability density of the time t of the next nucleation event is —dPo/dt. It follows that the mean value of 
Po at the time of nucleation satisfies 

P = l° O d tP (t)(-f)=l-P (16) 

and thus Po = 1/2. In the numerical implementation, we therefore monitor the increase of Po during the growth of 
the top terrace and create a new top terrace when Po = 1/2. For the nucleation rate uo the expression (H) will be used. 

C. APPROXIMATE ITERATIVE SOLUTION 

A full analytic solution of the model is difficult because the nucleation probability (|l5|) depends on the size of the 
top terrace, which is determined by the entire history of the wedding cake through the coupled equations (|l3|). Some 
first insight can be gained by assuming that the former top terrace ceases entirely to grow once a new island has 
nucleated on top of it. This should produce a lower bound on the island radii reached at a given time. 

Denote by P^°?i the radius of island n — 1 at nucleation of island n, and by t n the time of this nucleation event. 
Then during its tenure as the top terrace, the radius of island n grows, according to (|l4j), as 



Rn{t) = V^C* (17) 

and the time i n +i of the next nucleation event is determined by 

/ ir 2 F 2 (R top ) 5 r tn+1 \ 
Po(Wi) = exp TTT^- / dt [F(t - t n )f' 2 = 1/2. (18) 



2v> 



This implies a recursion relation 



R^^RTiR^f 7 (19) 
for the size of the top terrace at nucleation of the next layer, where we have introduced the characteristic radius 

71n2\ 1/5 /*/\ 1/5 (v>^ 



P C =^— j y-) , 0.868.^-j . (20) 

The time interval t„ = i n +i — t n during which nt op = n is related to P^ op through ([l7|), and correspondingly satisfies 

t^^F-^Ft^ 7 . (21) 

It is easy to check that the recursion relations ([l9|) and ( f2l| ) approach fixed point values 

P4 op -> Rc, Fr n -> 1 (22) 

exponentially fast in n. Thus asymptotically there is one nucleation event during the growth of one monolayer^, 
and nucleation occurs when the radius of the top terrace has reached the value R c . We shall see below that these 
statements remain valid for the full dynamics, although the approach of Pj; op and r„ to their asymptotic values is 
slower than exponential due to the coupling to the lower layers. 



3 This has also been observed experimentally for engineered mounds grown on Cu(100) 
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D. NUMERICAL INTEGRATION 



Since the deposition flux F and the size i?o of the base of the mound can be eliminated by rescaling the time and 
the island radii, the model effectively contains only the single parameter v' /F. For convenience we will express this 
parameter through the coverage 9 C of the top layer at nucleation, which is given by 

9c -{r~J fv^ - (23) 



, i 1 1 i , r 
i 1 i / 

i I 

i ~n 



j i 

i 1 ■■ , / A, 

i 1 ■■ , / -A, 

,■■ \ -A. 



FIG. 3. Wedding cakes grown numerically with 6 C = 0.1, at total coverages of Ft = 5, 10 and 15 monolayers. The figure 
shows a stack of circular islands seen from the side, i.e. the bar at level n has a length proportional to 

Figure 4 shows the result of a numerical integration of the model equations for 9 C = 0.1. The figure suggests that 
the number of exposed layers (with coverages < 9 n < 1) increases with time. This is quantified in Fig. 5, which 
shows the evolution of the surface width 



W 2 (t) = Y,{n-Ft) 2 {9 n ^x-e n ). 



(24) 



In the case of purely statistical growth with v' = 9 C = 0, W 2 = Ft because the exposed coverages n _i — 9 n are given 
by a Poisson distribution with parameter Ft p| , p6| . Figure 5 shows that also for 9 C > the squared width grows 
linearly in Ft, but with a coefficient that is less than unity; in Section III E| it will be shown that the coefficient is 
simply 1 — 9 C . Since the lateral mound size is fixed, the unbounded increase of W(t) implies that the mounds steepen 
indefinitely, with the terrace width on the mound slopes decreasing as l/\/Ft; a quantitative expression is given in 
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FIG. 4. Surface width ( pi| ) as a function of time. 
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Finally, Fig. 6 illustrates the convergence of the top terrace size at nucleation, and the time intervals between 
subsequent nucleation events, to their limiting values predicted by (^2|). As might be expected from the behavior of 
the surface width, the convergence is linear in 1/y/Ft, rather than exponential. 
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FIG. 5. Convergence of the time between nucleation events, and the top terrace size at nucleation, as a function of \jyfFt. 



E. THE ASYMPTOTIC MOUND SHAPE 
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FIG. 6. Convergence of numerically generated wedding cakes to the asymptotic shape function <E>(a;). The symbols show the 
coverage profile rescaled according to (Est) . 
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The numerical solution of the model equations shows that the height profile of the mound converges to a time- 
independent asymptotic shape, when viewed relative to the mean film height Ft and rescaled horizontally by y/~F~i 
to compensate the increase of the surface width (Figure 7). To derive the shape analytically, we write the coverage 
profile 9 n (t) as a traveling wave of width yFt, 

e n {t) = $[{n-t)/sJTt]. (25) 

Inserting this into (|]1]) and expanding for large Ft one finds that the shape function $(x) has to satisfy the differential 
equation 



= -x&{x), 



(26) 



where x = (n — t)/y/Ft is the scaling variable. This shows that the inflection point of the profile, where <£>" = 0, is 
always located at x = 0, i.e. at n = Ft. The solution of ( fzlf ) which satisfies the boundary condition lim^^-oo $(a;) = 1 
reads 



$(x) = 1-J-C dy e~ y /2 = 1 - C[l + eri(x/V2)], 



(27) 



where erf (s) denotes the error function, and C is a constant of integration. The profile is cut off at the rescaled height 
£max of the top terrace, where the coverage takes the value 9 C , 



$(lmax) = 9c- 

Accordingly, the height of the top terrace above the mean film thickness Ft is x max ^/F~i. 



(28) 
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FIG. 7. Scaled asymptotic mound shapes for 6 C — 0.01, 0.15, 0.4 and 0.8. 
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FIG. 8. Shape parameters C and x ma x as a function of 9 C . 

The two parameters C and x max of the shape function are related by (|28|) , but to fix both of them a further relation 
is required. This is provided by the normalization condition 



Ft, 



(29) 



which, using (|25|), translates into 



dx <&{x) 



dx (1 - ®(x)). 



(30) 



Together equations ( |27| , |28| , |30D define a family of shape functions parametrized by 8 C . Examples of these functions are 
shown in Fig. 8, while Fig. 9 illustrates the dependence of x max and C on 9 C . For the case of Pt growth on Pt(lll) 
at 440 K, the experimentally determined mean mound shape has been compared to the scaling function $(x) for 
8 C — 0, and good agreement was found except near the top of the mound |]Io| . Evaluation of the top terrace size and 
comparison with ( |20| ) yields the estimate AE's « 0.25 eV psj . 

Other features of the mound shape can be easily extracted. The surface width turns out to be given by 



W 



(1 - 9 c )Ft 



(31) 



asymptotically. The mound radius at the level of the mean film thickness, x = 0, is given by \J\ — CRo, and the 
maximum terrace width on the mound slopes, located also at the inflection point x = 0, decreases with increasing 
film thickness according to 



'max /Ro 



$'(0) 



c 



1 



2 v /¥(0) VFi ^2tt{1 - C) VFi 



(32) 



F. THE DISTRIBUTION OF TOP TERRACE SIZES 



In an image of a multilayer film subject to the mound instability one will typically see top terraces of a range of 
sizes, because they are caught in different stages of their evolution. In jl7| a simple model for the size distribution 
was developed, which decouples the dynamics of the top terrace from the rest of the mound by assuming that each 
top terrace grows on a template - the base terrace - of constant radius Rb- The nucleation probability ( |l5| ) is then 
evaluated with 







-RtopW = VFtR b , 



(33) 



and its derivative —dPo/dt is interpreted as the probability distribution for the time intervals between two subsequent 
nucleation events. The base terrace radius is fixed through the requirement that the mean time between nucleation 
events equals the monolayer deposition time, which implies 



■L 



F / dt P Q (t) = 1. 



(34) 



One obtains R b « 0.867 • (v'/F) 1 / 5 , which is virtually indistinguishable from the expression ( po|) derived from the 
condition Pq = 1/2. If one demands instead that nucleation takes place at the maximum of the probability density, 
i.e. at the time determined by d 2 Po/dt 2 — 0, a prefactor of 0.873 is obtained. All reasonable criteria give the same 
result, because the probability density —dPo/dt is sharply peaked around its maximum. 
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FIG. 9. Distribution of sizes of top terraces, measured in terms of their perimeter L in units of lattice sites |L7|| . The full line 
shows a histogram obtained from the evaluation of STM images of 145 mounds, grown at 440 K in the presence of a partial CO 
pressure of 1.9 x 10 -9 mbar, at a total coverage of 37.1 monolayers. The thick dotted line shows the (appropriately normalized) 
prediction j35|). The thin dashed line shows the corresponding histogram for the base terraces, i.e. the terraces supporting the 
top terraces. It has been included to demonstrate that the assumption of a constant size of base terraces is reasonable. 

Within this model, the probability P(R) to observe a top terrace of size R is equal to the probability that the time 
since the last nucleation event is at least F~ 1 (R/Rf ) ) 2 , which equals Po(F~ 1 (R/ Rb) 2 ) ■ Changing variables from t to 
R then yields 



n 2 FR 7 



(35) 



As can be seen in Fig. 10, this gives a quite satisfactory description of the experimental data. By fitting the mean of 
(|35|), the step edge barrier for Pt(lll) in the presence of CO has been extracted from images of multilayer mounds [fL7| . 
The result AEs ~ 0.33 eV is in good agreement with the estimate of 0.36 eV obtained from second layer nucleation 
in the submonolayer regime (see Section [n|) . 



G. THE TOP TERRACE SIZE FOR REVERSIBLE AGGREGATION 



The derivation of (|20j ) can be repeated using the expressions (|9[|l0]) for the nucleation rate in the case of reversible 
aggregation, i* > 1 ]19[. In general the top terrace size is given by an expression of the form 



R c 



~F , 



(36) 
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where the exponents 7' and fi' take the values 



1 = = - 2 (regime II) 



(37) 



7 = 7i , o , A* = .« , o (regime III) 



i* + 3' r i*+3 



(38) 



In terms of the kinetic rates 1/ and F, the two regimes are defined by the inequalities Fjv <C v'/u -C {F/v) Sl 
(regime II) and {F/v) &1 < < {F/v) x/2 (regime III), where 



i*(2z* - 1) 
2i*(i* + 1) + 2 



and 



(39) 



X : 



z* + 2 



(40) 



is the island density exponent of classical nucleation theory which relates the density N of first layer islands to the 
ratio Fjv through p9] 



N 



(41) 



IV. CONTINUUM THEORY 



The continuum theory of mound formation is based on the notion of a growth-induced, uphill surface current 
P, p4||3f]j3l[ | . Figure 11 illustrates the basic idea: Because atoms deposited onto a vicinal terrace on the side of a 
mound attach preferentially to the ascending step, they travel on average in the direction of the uphill slope^j between 
the point of deposition and the point of incorporation. 

Fl 




FIG. 10. Illustration of the uphill surface current generated by step edge barriers. 

The evaluation of the current is particularly simple when interlayer transport is completely suppressed, and nu- 
cleation on the vicinal terrace can be neglected]^. An atom deposited onto a vicinal terrace of width I then travels 
an average distance 1/2 to the ascending step, and hence the current is Fl/2. Describing the surface profile by a 
continuous height function h(r,t) which measures the film thickness (in units of the monolayer thickness) above a 
substrate point r, the local terrace width is I = |V/i| _1 , and hence the current is given by the expression 



J = 



F 



rV/l. 



2\\7h\ 2 

The surface profile evolves according to the continuity equation 



(42) 



4 It is important to realize that this does not require atoms to actually climb uphill across steps 
5 Expressions for more general situations are derived in Ej|Jjl] |. 
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f + V-JW. (43) 



Here we are specifically interested in radially symmetric mounds. Rewriting ( [43|) in polar coordinates and using ( |42| ) 
yields the following evolution equation for the mound profile h(r.t): 



dh F d fdh\ 1 



dt 2r dr \ dr 



F. (44) 



We want to show that fl44[ ) possesses solutions corresponding to the asymptotic mound shape derived in Section |IIIE 
To this end we make the ansatz 



h(r, t) = VFt ip(r) + Ft. (45) 
Inserting this into (|44|) yields the differential equation 

rip(r). (46) 



dr \tl>'(r) 

In order t o esta blish the equivalence between the shapes described by ip{ r ) an d the scaled coverage distribution <I>(a;) 



of Section HIE , we need to verify that the function tp(r) defined implicitly by 

2 



r 

Rq 



$0(r)) (47) 



solves (|4^). Indeed, taking the derivative with respect to tp on both sides yields 



dr J 2 

which reduces to ( |46| ) upon taking another derivative with respect to r. 

This calculation shows explicitly that the sides of the mounds evolve according to the continuum equation 
in |ll| the same result was obtained for a one-dimensional geometry. The continuum description does however not 
include the nucleation on the top terraces, which have to be added to the profile determined by (^) as a boundary 
condition. The development of a continuum theory of epitaxial growth which explicitly incorporates nucleation remains 
a challenge for the future. 
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